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1. Introduction. 

Massless integrable quantum field theories in 1+1 dimensions and the associated fac- 
torized "massless scattering" have been intensively studied recently [jl[, ||. The general 
context where such theories appear is when the perturbation of a conformal field theory 
by a relevant operator causes it to flow to another conformal field theory. In the midst of 
the flow, the model has massless excitations but is not scale-invariant. 

In ref. || it was shown that the sine-Gordon model with imaginary potential (and its 
natural generalizations to higher level and rank) plays a fundamental role in understand- 
ing flows between minimal conformal field theories, like the usual sine-Gordon model for 
massive perturbations. Exact (LL) and (LR) S matrices were found phenomenologically 
for this problem using the Yang Baxter equation, unitarity and crossing-symmetry . The 
physical content of the model and the massless flows were then analyzed by thermodynamic 
Bethe ansatz (TBA). 

The purpose of this letter is to derive these S matrices from the exact solution of the 
quantum field theory. It is easy to check that the known equivalence between sine-Gordon 
and Thirring extends to our case. We shall therefore study the massive Thirring model 
(MT) with imaginary mass (IMT), which is easier to handle technically. We use the well 
known coordinate Bethe ansatz (BA) 0] , build the physical vacuum of the model and show 
explicitly that there is no gap in the excitation spectrum. The model possesses non-trivial 
(LL) and (LR) scattering. The corresponding S-matrices are found from the BA equations 
and shown to be in agreement with the results of . 



2. Structure of the model. 

The massive Thirring (MT) model, defined by the Hamiltonian 

H = J dx[-i(^+d x ^ - ip}d x ifo) + mo(^ + V^i) + 2^+V>+V 2 M (2-1) 

is a well-known example of an exactly soluble model. Its quantization on a circle of length 
L, first carried out in Q, leads to the Bethe ansatz equations for the rapidities of N 
pseudoparticles: 
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Here p(0) is the momentum of a particle, — 0j) is the two-particle bare phase shift, 

, sinh±(2zu- 0) 1 , . 

and Ii is a set of iV integer numbers. The momentum and energy h(/3) of a particle lie on 
the hyperbola h 2 — p 2 = m 2 ,, which is parametrized in terms of rapidities a: 

p = mo sinh a, h = mo cosh a. 

Suppose now that the mass becomes purely imaginary, mo = im. The hyperbola gets 
rotated, p 2 — h 2 = m 2 , and momentum exchanges its role with energy. For the energy 
and momentum to be real, the rapidity should now be a = — + 0, where belongs to 
{Im0 = 0} |J{Im,9 = n}. The new parametrization in terms of 3 is 

p = m cosh 0, h = m sinh 0. (2-4) 



This affects the solutions of equations ( |2.2|) and modifies drastically the structure of the 
ground state and excitations of the Thirring model. Note that the presence of an imaginary 
mass in the Hamiltonian ( |2.1| ) makes it non-hermitian. However it is probable that the 
spectrum is real. Indeed complex values of are strongly restricted by the Bethe equations. 
The latter allow only string type solutions with a definite pattern of complex 0's, for which 
the total energy is real. 

In the thermodynamic limit, equations ( |2.2| ) become coupled integral equations for 
the set of densities Pj{0)- Let us discuss first the ground state of the model. The densities 
of the ground state distribution p± satisfy the following system of integral equations (dots 
denote derivative with respect to 0): 

/B />oo 
<&i,i(/3 - a)p+(a)da - / $- lfl (0 - a)p-(a)da - 2irp + {0) 

7 (2-5) 

/is poo x * 

$_i 5 i(/? - a)p + (a)da - / - a)p-(a)da - 27rp_(/5) 

-oo J-B 



where + iir) = = (see (|Q|)). The subscript + (resp. -) means 

that the corresponding function is defined on the axis Im/3 = (resp. Im/3 = n). B is a 
cut-off such that p + = for > B and p_ = for < —B. Equations ( |2.5| ) define a 



one-parameter family of solutions p±(fl,B). The ground state will be found in principle 
by requiring that p± > and that the total energy, 



E = / + / h+p+, (2.6) 

J-B J-oo 

is minimal. One also has to check that this ground state is stable, ie that all excitations 
have positive energy. Note that the presence of pseudoparticles with positive bare energy 
in the ground state cannot be ruled out at this stage. The filling of such states brings some 
positive contribution to £0, but it affects also the distribution of all the other pseudopar- 
ticles through the Bethe equations (|2.2| ), an effect which can lower the total energy. 

Let us comment on the sign of the last term at the right hand side of ( |2.5| ). This 
sign depends on whether the integers Ii in ( |2.2| ) are an increasing or decreasing function 
of the bare rapidity. In other words it indicates the signature of the derivative y((3) 
where p(f3)L + ^ ■ — (3j) = y((3) (assuming this is /3-independent) . In the simplest 
situation, when the signs of p((3) and J2j ~ Pj) coincide for all (3 one can make a 
rigorous conclusion about the signature of y. This was the case in the usual Thirring 
model. Otherwise, the evaluation of sign(y) becomes more delicate because of the infinite 
summation over the so far unknown positions f3j, the parameters L, N and the momentum 
cutoff which should be sent to infinity. Close to the free point one expects however that 
the sign of y is the sign of p for the rapidities of filled states; this is how we choose the 



sign in equations (|2.5| ) . Our choice gives a reasonable solution of ( |2.5| ) for which p± > 
and B = 00. Other choices seem to lead to negative densities, but it is difficult to argue 
for finite B except close to the free point. 

It is known that in order to solve the MT model appropriate regularization is required 
0] . To eliminate the divergences we introduce an analytic momentum cutoff similar to the 
one used in the £/(l)-symmetric Thirring model with two isospin degrees of freedom || 
and in the XXZ lattice model ||, 0]: 

. . cosh|(^- f-ip) . sinhi(/3 + / - ip) 

p+{p) — im log f- + im log f- - 

cosh 2 [P — f + ip) sinh ^ [(3 + f + ip) 

= G((3-f) + G((3 + f + i7r), (2.7) 

1 1 



p+((3) /'/sin// 

\cosh(p — j) + cos p cosh(p + j) — cos p t 

and P-(f3) = p+((3 + in). The peculiarity of the cutoff ( |2.7|) is an unusual shift of ar- 
gument by in between the first and second terms. Upon the cutoff removal, / — > 00, 
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we recover our momentum of pseudoparticles: p + — > Am sin/xexp(— /) cosh ,9 , p_ — > 
—Am sin p exp(— /) cosh (3. 

In principle, the system of integral equations (|2.5| ) cannot be solved analytically for 
arbitrary B. Using the symmetry p+(j3) = p-(—(3), that follows from p_(—f3) = p + (f3), 
the system ( |2.5| ) can be rewritten as a single Fredholm integral equation: 

MP) = ~ I - «) + ^-i,i(/3 + a)) p+(a)cfa - 2tt P+ (/3). (2.8) 

Equations of such a form have previously appeared in related problems, for instance the 
£/(l)-symmetric Thirring model in a magnetic field (WC sector) ||, @. An approximate 
solution can be obtained if one treats the second kernel in (|2.8| ) as a small perturbation to 
the first one and then makes use of the Wiener-Hopf method. 

We show that in the repulsive regime (0 < p < 7r/2) the solution of (|2.5|) with B = oo 
is consistent, and we conjecture that it represents the ground state of the model at least 
for < p < n/A (strong repulsive regime). Our main argument is that such a solution 
leads to reasonable density distributions, to positive excitation energies over the physical 
vacuum, and to the correct scattering matrices (see section 3). Solving equations ( |2~5| ) 
with B = oo by Fourier transformation and using 

sinh(7r - 2p)k f°° ■ ( ^AkB-, a sinh2/i/c 



f $i 1 ((3)e^d(3 = 27r Sinh(7T ~, f *-i i{PV kp d(3 = -2tt 

J-oo ' smhyrA; J_ 00 

P ^e^d(3 = me*^ - me^^'^. 
27T sinhnk sinh7r/c 

we obtain: 

'p + (k)\ m ( e~ lk f 



sinh7r/c 



y p-{k)J 2 cosh pk \ e lk f 
The inverse Fourier transform p(f3) = J p(k) exp(—ik(3)dk/2% yields: 



(2.9) 



m 171 
P+ = 4^cosh^(/3 + /)' P ~ = 4/.cosh^(/3-/)- ( ' } 



Upon cutoff removal, 



m -tP* 
P±\P) ~^ TT e 2 »e +2 », f — > oo. 



2(jl 

The density p + (p_) decrease exponentially for those (3, where h{(3) 3> m, i.e. (3 — > +oo 
(f3 — > — oo). On the other side of the deep ultraviolet region, (3 — ► — oo ((3 — > +oo), where 
h(i3) <C — m, the functions ( |2.10| ) have the same asymptotic behaviour, as the density in 
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the usual MT model ||. This is expected since for (3 — > — oo (f3 — > +00) we can neglect 
the contribution from the pseudoparticles located on Im(3 = 11 (Im/3 = 0) line in eqs. ([2.2p . 
When fx — > 0, P±(P) ~ Hfafi), where i?(a;) is the Heaviside step function. 

To obtain the energy spectrum of excitations at T = we use the TBA method. 
Variation of densities of pseudoparticles over the physical vacuum leads to the system of 
integral equations ||: 

hj = elf + A jk * Ok£~k 

J _ (2.11) 

Pj = tt/ + Ajk * a k n k , 

where * denotes convolution, and 7r+ are the energies and momenta of different excita- 
tions labelled by j (e + > and e~ < by definition), and 

A jk = -$jfc/27r + (TjdjkSifi), 

where aj = — 1 for the holes (j = ±1) and aj = 1 for \j\ > 2. For technical simplicity we 
have restricted the coupling parameter to rational values \i = 7t/s, s being an integer. In 
this case the allowed strings are those of length from 2 to s — 1, centered around the axis 
Im/3 = 0, 7r which we denote by ±j (strings and antistrings in the usual denomination). 
To see this one follows the usual argument, taking the modulo of the both sides of ( |2.2j ) 
0. When there is a real l3\ in the left hand side of (pT2|), the modulo of both sides is 
easily seen to be 1. When the rapidity of one of the particles of the string is plugged into 
the L.H.S., one of the terms in the R.H.S. becomes equal to zero or infinity. The L.H.S. 
vanishes or blows up as well in the limit L — ► 00, depending on the sign of the imaginary 
part of the momentum. The latter is however a cumbersome expression depending on the 
cut-off /; this is not usually the case. If one takes the limit / — > 00 at first, one reproduces 
the expression for the momentum ( |2.4|) , from which the result follows straighforwardly. 

We choose the regularization of the bare energy to be hj (l3) = pj (l3) , see ( |2.7|) (another 
regularization, such as hj([3) = G(j3 — f) — G((3 + f + itv) could lead to different mass 
renormalization, which in our case is not important). 

For the holes in the ground state (|2.10|) , solution of (|2.11| ) with e^j = yields in the 
limit / — » 00: 



These excitations are massless left and right-moving particles which we call by analogy 
" solitons" . 

For the strings of length \j\ = 2, s — 2 we obtain £ j~ = 7T j~ = 0, where we used 

- cosh aksinh( tv — j u)k ~ .,. coshn/csinh ink 

simvnk sinh7r/c 

For the longest (s — 1) string the calculation yields ^(s-i) = k^il- Such a string with the 
center located on the line Im/3 = (Im/3 = n) behaves as a hole with the same rapidity 
located on the line lm.fi = n (Im/3 = 0). 

There is one more peculiar excitation in the IMT model in thermodynamic limit. 
This is a 1-string, located on the imaginary rapidity axis exactly in between the two filled 
vacuum lines (see fig. 2), with (3 = in/2. To check that j3 = in/2 is indeed a solution of 



( | . 2| ) in the thermodynamic limit we take the modulo of both sides of (|2.2|) , where the 
momentum is given by (|2.7|) , and set L, f,N — > oo. If one sets /3j = x + in/2 in the left 
hand side and combines the terms with f3j and /3j> = 2x — f3j + in in the R.H.S. of ( |2.2| ), 
the product of terms in the R.H.S. gives 1. This can be seen also by symmetry of the 
picture, depicted in fig. 2. In the L.H.S., the momentum in the exponential is proportional 
to mexp(— /). In order to keep the physical (dressed) mass finite, one should introduce 
the cut-off dependence in the bare mass, m ~ exp(^). Then, in the limit L — > oo the 
modulo of the L.H.S. is 1 iff x = 0, independently from /. If one chooses f3i real in L.H.S. 



of ( 2.2|) , then f3j = in/2 will appear in the R.H.S., making its modulus to be different 
from 1. So it seems f3 = in/ 2 is not allowed. However, observe that multiplying this 
equation by ( |2.2|) with j3i> = —fli + in we obtain 1 in the R.H.S., too. This suggests that 
the Bethe equations describing such 1-string should be combined into pairs according to 
the symmetry shown on fig. 2. The possible physical explanation is that the vacuum of the 
IMT model consists of a condensate of pairs of L and R moving particles with opposite 
momenta and equal energy. The appropriate Bethe equations would then represent the 
periodic boundary conditions, written for the wave function when one moves such a pair 
through the system. The formation of pairs appears actually necessary to conserve parity, 
which is not broken in our problem. 

The physical (dressed) energy and momentum of this new excitation, calculated by 
means of ( |2.11|) vanish. 

'0 — ^0 



where we used the phase shifts ( |2.3|) given by 

\<P_1 0/ smhnk \ e 2 J 

and 

, ... , smh^k far , .... 
^ (fc)=4 "^mh^ C ° Sh( T + Zfc/) - 

This result still holds formally for excitations /3 = a: + if, x real. 

We conclude this section by remarks on the structure of the IMT model in the attrac- 
tive regime and in the vicinity of the free fermion point in the repulsive regime. When one 
goes to the attractive regime (fx > tt/2), the two-particle bare phase shifts ( |2.3|) change 
sign: — > —<&(/?) when |U — > tt — /i. The solution of ( |2.5[ ) with S = 00 is no longer valid 
since it leads to negative densities. Therefore, B must be finite. It is natural to suppose 
that a phase transition happens in the system at \x = tc/2 or at \x = A < tt/2, where 
A takes one of the possible values n/3 or n/4. Such a phase transition would probably 
involve level crossing, like for the MT model in the repulsive regime and appropriate cut 
off, see §. 

In the limit (jl — > n/2 there are two different solutions to (|2.5|) , obtained for B = 00 
and B = 0. The first one is given by ( |2.10|) , and the second is 2irp± = —p± (defined on 
the appropriate half-line). The ground state energies ( |2.6|) for these two solutions turn out 
to be identical, which suggests that the phase transition takes place exactly at \i = tt/2, 
and that the free point is singular, as anticipated in 0. A rigorous mathematical analysis 
of the integral equation ( |2.8|) or a detailed numerical solution (which is delicate because of 
cut off effects) could shed some light on this question. 

As in the real mass case § we expect that the physical results depend on the type of 
cut-off. It is possible that a very different behaviour may emerge for instance if one uses 
a sharp cut-off such as p±(j3) = if \(3\ > A. 



3. Scattering of excitations. 

In order to compute elastic S-matrices we adopt the method developed in 0. Let us 
consider holes at first. By definition, the phase shift for scattering of two holes is given by 

S h (p 1 ~(3 2 ) = -\ogS = ^ 1 -^ (3.1) 

% 



where <pi is the phase gained by a hole when going around the system and p 2 the same 
phase but in the presence of another hole. The are composed of a sum of two-particle 
bare phase shifts; for example 

ip 2 = L P± {p x ) + - Pi) + E - /%)• 

Here the sum is taken over all the solutions /3j of the Bethe equations (|2.2|) with two 
holes at positions j3\ and fa. The result of subtraction (pi — ip 2 is proportional to the 
backflow function of vacuum: 5h = 27rF(fa\fa). The latter is defined using the difference 
of the two solutions of the Bethe equations obtained with and without a hole at fa, 
F{fa\0) = (P ~ (3)Lp(f3). In our case we have two backflow functions F + , F_, describing 
a backflow on Im/3 = and Im/3 = tc lines respectively. One can show using ( |2.2|) that F± 
satisfy the following system of integral equations: 

- fa) = Ki * F + + $-1,1 * + 2ttF + 

(3-2) 

$_i,i(/3 - A)) = $1,1 * F- + $-1,1 * ^+ + 2ttF_ 

Equations ( |3.2|) describe the backflow caused by a hole at j3 = fa on the real axis; for a 
hole on Im/3 = n axis simply substitute F + <-> F_ in eqs ( |3.2| ). Taking a derivative with 
respect to (3 and applying the Fourier transform to both sides of ( |3.2| ) we arrive at the 
following solution (in Fourier space): 

F+\ _ e lk/3 ° fsmh(n -3fi)k\ 

(3.3) 



Ft/ 2 cosh/i/c sinh(7r — 2fj)k \ — sinh^/c /' 

i*+ represents scattering of two holes on the same line ( ie LL or RR scattering) and F- 
scattering of two holes on different lines ( ie LR or RL scattering). From ( |3.1|) we obtain: 

- W) . (3,, 

i ^ ogAMe) = - /;; ^ 2cosh ^_ 2rt ^ , (3.5, 

where 9 = fa — fa. Note that this method does not fix the constant normalization factor in 
the matrix elements All-, Arl, which should be fixed by other constraints (e.g. unitarity). 
Eventually, we obtain: 

A (9) * [ + °° ' k9 s ^ n ^ 1 ( 7r ~~ 3//)/c cZfc ^ 

2 J-oo cosh/x/c sinh(7r — 2^u)/c fc ' 



... i f +OG , „ sinhn/c <i/c , „. 



Expressions ( |3.6| ) and ( |3.7| ) are in agreement with results of || , which enables us to identify 
coupling constants in imaginary mass MT model with those of imaginary potential SG 
model: /u = ir — /3|q/8. Note that the integrand in (|3.6| ) and ( |3.7|) blows up when fx — ► 7r/2 
(free fermions). The LL-scattering becomes trivial when — > 7r/3, whereas for the LR- 
scattering this happens when fi — > 0. The L and R systems are decoupled at = and 
the model not only has massless excitation spectrum, but also is scale invariant. 

In the same manner one can calculate the backflow functions for strings. We will need 
them to get the rest of the matrix elements for the soliton-antisoliton S-matrix. For the 
j-string (centered around the real axis) one has a system, similar to (|3.2|): 



(3.f 



- $ ld (p - f%) = $1,1 * F+ + $_i,i * F_ + 2tvF+ 

- $-ijCS - A)) = $1,1 * F- + $-i,i * F+ + 2ttF_ 

The solution to these equations for j < s — 1 is: 

ps = sinh(7r- (j + l)n)k clk0a p s = sinhjj -l)fxk cik0a g ^ 
+ sinh(7r — 2fi)k ~ sinh(7r — 2fi)k 

and for j = s — 1 

ps = sinh^fc clk0o ^ ^ s = sinh(7r-3/x)fc 

+ 2cosh / u/csinh(7r — 2fx)k 2cosh / u/csinh(7r — 2/j)k 

which is the same as the solution for the hole on Im/3 = 7r axis. 

To address the question of soliton antisoliton scattering it is useful to determine the 
charge of excitations, i.e. the eigenvalues of the operator N = J ipiipi + V^V^- Setting 
the charge of the vacuum to zero, one has 



/+oo 
(Fj + F^)d[3 
-oo 



2(tt-2/z) 



(3.11) 



Q str =j+l (Fl + F s _)dj3 = ' (3.12) 



oc 



for a hole and 



for the j-string, where j < s — 1 and ( |3.3|) , ( |3.9|) were used. The simplest configuration of 
vanishing charge can be composed of two holes and a 2-string (see fig.l). It reproduces LL 
(RR) soliton-antisoliton S-matrix. Indeed, one can see using the definition ( |3.1] ) that the 
additional phase which the hole {j3\) acquires in the presence of the string + /?2)/2) 
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and another hole (fa) consists of the sum of two phase shifts of a hole- hole and string- 
hole scattering: dhsh(0) = $h(0) +£fc«(0/2). The hole-hole scattering phase was calculated 
previously (see (|3.6| )). As for the string-hole phase shift, it is given by the solution of 
equations (|3.8| ) for j = 2 and fa = +/32) /2, F + and F_ being related to the transmission 
and reflection amplitudes. Taking the inverse Fourier transform of (|3.9| ) we obtain after 
some manipulations: 



sinh^^ _ cosh |( » 



As in the usual Thirring model, solutions of the Bethe ansatz equations correspond to 
linear combinations of multiparticle states for which the scattering is diagonal. In the 
usual basis the matrix elements would be All (describes ss — > ss process), Bll (describes 
sa — > sa process), and Cll (describes sa — > as process), where Bll ± Cll = AllU^ l . 

To describe the LR (RL) scattering we introduce two holes and a 1-string at (fa + fa + 
z7r)/2, formally without the constraint that the sum of rapidities vanishes (as illustrated 
on fig. 2 for fa = —fa). The backflow due to our new 1-string can be obtained from the 
system (|3.8|) . The result is 

( P +\_ JKth+M/2 sinh/xfc (eW>\ 
The charge of this 1-string can be evaluated by formula ( |3.12| ) to be: 

7T 

Qiir/2 = (3-15) 

TV — l\l 

so that the total charge of our state is Qi w /2 + 2Qh = 0. As in the case of LL-scattering, 
the functions ( |3.14[ ) lead to the string-hole phase shifts which are related to the LR 
transmission-reflection amplitudes. We get: 



e\ n k^(—Q- v Jilt i\ rnshlLf— i I — j) 

blIln 4 U-2 M J_ 7i— 2 M V j J+ __ CObX1 4 U-2/x J_ tt-2h L ) 

sinh^f— e - ^--{\ RL ~ cosh^f— e - - i) 

° um 4 W-2/Lt tt-2 m ^ LUbI1 4U-2/x 7r-2/Lt *7 



The hole-hole scattering is given by ( |3.7| ). Putting everything together, we obtain the LR 
soliton-antisoliton S- matrix elements Arl, B^l and Crl, where Brl ± Crl = ^-RlU^ l . 
This S-matrix is identical to the one found in by solving the YBE. One should remember 
however the constraint on the position of the 1-string, obtained at the end of section 2, that 
imposes to consider only the scattering of pairs of L and R moving particles with opposite 
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momenta. Moreover note that due to the presence of purely exponential factors in the 
Fourier images F±(k), the functions F±(f3) have imaginary parts (however, the combination 
F+ + F- is a real function). This means that particles of the vacuum are pushed off 
the Im.p = and Im/3 = 7r lines in the complex rapidity plane towards the particle at 
+ @2 + nr)/2. So it may be that the complete scattering theory is more complicated 
than the one of || . Probably a numerical study of the appropriate inhomogeneous 6- vertex 
model could shed light on this. 

Note that there are misprints in the formulas (2.7) of [[|, page 5757. The correct form 
of these equations is given in the preprint version of . 
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Figure captions 

Figure 1: a pair of holes and a 2-string centered on the real axis describe a pair of L moving 
" soliton" and " antisoliton" . 

Figure 2: a pair of holes with opposite rapidities on the real axis and the axis lm/3 = n 
together with a string at (3 = in /2 describe a pair of L and R moving "soliton" and 
" antisoliton" . 
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This figure "figl-l.png" is available in "png" format from: 
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This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9403022vl 




